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ABSTRACT 

Wc derive inner dark matter halo density profiles for a sample of 165 low-mass galaxies using rotation 
curves obtained from high-quality, long-slit optical spectra assuming minimal disks and spherical 
symmetry. For p(r) ~ r~ a near the galaxy center we measure median inner slopes ranging from 
a m = 0.22 ± 0.08 to 0.28 ± 0.06 for various subsamples of the data. This is similar to values found by 
other authors and in stark contrast to the intrinsic cusps {a.i n t ~ 1) predicted by simulations of halo 
assembly in cold dark matter (CDM) cosmologies. To elucidate the relationship between a m and a^t 
in our data, we simulate long-slit observations of model galaxies with halo shapes broadly consistent 
with the CDM paradigm. Simulations with cti n t = 1/2 and 1 recover both the observed distribution 
of a m and correlations between a m and primary observational parameters such as distance and disk 
inclination, whereas those with oti n t = 5/4 are marginally consistent with the data. Conversely, the 
hypothesis that low-mass galaxies have a^t = 3/2 is rejected. While the simulations do not imply 
that the data favor intrinsic cusps over cores, they demonstrate that the discrepancy between a m and 
ai n t ~ 1 for our sample does not necessarily imply a genuine conflict between our results and CDM 
predictions: rather, the apparent cusp/core problem may be reconciled by considering the impact of 
observing and data processing techniques on rotation curves derived from long-slit spectra. 
Subject headings: dark matter — galaxies: dwarf — galaxies: halos — galaxies: kinematics and 
dynamics — galaxies: structure 



1. INTRODUCTION: LONG-SLIT SPECTRA AND THE 
CUSP/CORE PROBLEM 

Within the standard cosmological cold dark matter 
(CDM) paradigm, structures form hierarchically via 
gravitational collapse of primordial density fluctuations. 
This framework has been successful at providing physi- 
cal interpretations for a wide range of phenomena, no- 
tably the angular power spectrum of cosmic microwave 
background anisotropies, the large-scale structure of the 
galaxy distribution, and fundamental scaling relations in 
disk galaxies. 

In recent years, significant theoretical progress has also 
been made in predicting halo shapes obtained from sim- 
ulations of structure formation in a (A)CDM universe. 
Detailed collisionless numerical simulations of halo as- 
sembly (see Reed et al. 2003 and references therein) have 
led to a general consensus on the fundamental properties 
of the resulting CDM halos: their density distributions 
p(r) deviate significantly from the power laws predicted 
by analytic calculations (e.g. Gunn & Gott 1972); they 
have a broadly "universal" shape over many decades in 
mass that can be parameterized by simple fitting formu- 
lae; and they are cuspy, in that p(r) cx r~ aint at the 
smallest halo radii r probed by simulations, and the in- 
trinsic inner slope aint is of order 1. 

A simple analytic profile that encompasses these three 
characteristics was first proposed by Navarro et al. (1996; 
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1997, hereafter NFW) for dark matter halos: 

= i i / M2 ' W 

{r/r s )[l + [r/r s )Y 

where p s and r s are the characteristic halo density and 
radius, respectively. Subsequent generations of simula- 
tions showed some disagreement over the value of tti n t 
that best described the resulting halos: inner slopes shal- 
lower than the NFW value of a^t = 1 (Subramanian 
et al. 2000; Taylor & Navarro 2001; Ricotti 2003) and 
as steep as aj„t=1.5 (e.g. Fukushige & Makino 1997, 
2001; Moore et al. 1998, 1999) have been proposed. A 
systematic study of these discrepancies has been under- 
taken by Navarro and collaborators, using a suite of high- 
resolution simulations that reliably probe galaxy-size ha- 
los down to small fractions of their virial radii (Hayashi et 
al. 2004; Navarro et al. 2003; Power et al. 2003). They 
find that while the density distribution does not converge 
to a well-defined asymptotic power-law at the innermost 
resolved radius, the average a^t is 1-1, 1.2, and 1.35 
at 1% of the virial radii of cluster, galaxy, and dwarf 
halos, respectively. The NFW slope a^t = 1 remains 
consistent with the simulated halo shapes, but cusps as 
steep as a^t = 1-5 are ruled out in almost all cases. 
Some impacts of the presence of baryons on these halo 
shapes have been investigated (e.g. Mo & Mao 2004); 
however, fully self-consistent simulations of halo shapes 
including a variety of gasdynamical effects have not yet 
been achieved. 

While collisionless halo shapes are consistent with the 
average properties of high redshift halos inferred from 
weak lensing studies (e.g. Hoekstra et al. 2004), on 
galaxy and cluster scales the agreement between theory 
and observations is much less certain. The majority of 
observations rely on luminous tracers embedded at small 
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Fig. 1. — Sky distribution of sample galaxies. Solid red lines 
denote the plane of the Milky Way (b = 0° ± 20°). 



r to probe the halo potential. Even though some stud- 
ies of halo shapes prefer CDM profiles (e.g. Mahdavi & 
Geller 2004; Pointecouteau et al. 2004), a variety of evi- 
dence suggests that at least some real halos do not have 
the cusps predicted for the pre-baryon collapse systems, 
but rather cores in which the inner halo density has only 
a weak dependence on r (or none at all). For instance, 
studies of gravitational arcs around clusters (Sand et al. 
2004; but see El-Zant et al. 2004), models of gas flow 
along strong bars (Weiner et al. 2001), and the mass 
budget in the Milky Way from microlensing and stellar 
kinematics (Binney & Evans 2001; Binney 2004) all yield 
measured inner halo slopes a m that are much more "core- 
like" (a m < 0.5) than the cusps predicted by CDM. The 
rotation curves (RCs) of spiral galaxies have also been 
extensively used to constrain halo structure. In particu- 
lar, Salucci (2001) concludes that halos have cores larger 
than the corresponding disk scale-lengths from a robust 
analysis of 137 disk-dominated systems that is immune 
to biases on the inner RC shape. This apparent con- 
flict between CDM theory and observations in these and 
other systems has been called the "cusp/core problem". 

Among experiments that probe the structure of galac- 
tic halos, substantial resources have been devoted to de- 
riving inner dark matter halo shapes for dwarf and low 
surface brightness (LSB) galaxies. Although in some sys- 
tems the RC features correlate with those seen in the 
light distribution (Sancisi 2004), the global RC dynam- 
ics in dwarf and LSB galaxies are thought to be dark 
matter dominated even at small r (e.g. Swaters 1999). 
The advantage in inferring their halo properties from ob- 
served gas dynamics is therefore twofold: first, there may 
be little contribution to the observed RC from baryonic 
components, in which case the RC shape is a direct man- 
ifestation of the inner halo structure. Second, baryonic 
collapse in the halo gravipotential well during disk forma- 
tion may not significantly alter the halo structure. The 
halos of dwarf and LSB galaxies may therefore represent 
pristine, pre-baryon collapse distributions better suited 
to direct comparisons with collisionless CDM simulations 
than the halos of more massive or luminous systems. 

Early observations of dwarf and LSB kinematics from 
H I synthesis observations (e.g. Carignan & Beaulieu 



1989; Cote et al. 1991; de Blok & McGaugh 1997) yielded 
RCs that are well-described by models with cores rather 
than cusps. This apparent crisis for CDM (e.g. Flores & 
Primack 1994; Moore 1994) has since been called into 
question, however, because beam smearing caused by 
limited spatial resolution in the inner halo systematically 
lowers the measured inner profile slopes (Blais-Ouellette 
et al. 1999, 2004; Swaters 1999; Swaters et al. 2000; van 
den Bosch et al. 2000; van den Bosch & Swaters 2001; 
but see Gentile et al. 2004). Alternatively, long-slit opti- 
cal spectroscopy yields arcsecond-resolution RCs in rela- 
tively short integration times. There has therefore been 
a concerted effort to obtain long-slit spectra for many 
dwarf and LSB systems (e.g. Borriello & Salucci 2001; de 
Blok et al. 2001, hereafter dBOl; McGaugh et al. 2001; 
de Blok & Bosma 2002; Marchesini et al. 2002; Swaters 
et al. 2003, hereafter S03). A full characterization of the 
kinematics is then achieved by combining the resulting 
high resolution RCs with measurements of the outer RC 
shape in H I (e.g. de Blok et al. 2001; Marchesini et 
al. 2002; S03; Gentile et al. 2004). Approximately 70 
high-quality RCs have been obtained from these studies, 
and almost all are better described by cores than the 
cusps predicted by CDM. These data are clearly consis- 
tent with halos having cores rather than cusps; nonethe- 
less, mass models of hybrid Ha+H I RCs often cannot 
rule out the a int = 1 case (S03, Hayashi et al. 2004). 

Many studies have focused on a simpler measure of a m 
from long-slit spectra alone, in which the stellar mass is 
ignored and spherical symmetry is assumed (dBOl; de 
Blok et al. 2003, hereafter dB03; S03). This is clearly an 
over-simplification of the problem (e.g. Fall & Efstathiou 
1980; Salucci 2001), but one that produces meaningful 
results: minimum disks yield upper limits on the true 
steepness of the halo profile, and there is typically lit- 
tle change in the halo slope when a luminous component 
is introduced (dBOl, de Blok & McGaugh 1997). This 
technique, too, has revealed mostly core-like a m . How- 
ever, even at arcsecond resolution a variety of observa- 
tional uncertainties such as slit offsets, slit width and 
seeing, and geometric effects such as galaxy inclination 
all lead to an underestimate of a m relative to the un- 
derlying otint- Moreover, inner halo shapes are measured 
in the region where RCs are most susceptible to fold- 
ing errors, and where the relative uncertainties on the 
RC points are largest. As a result, raw long-slit spec- 
tra are heavily processed before the inner halo shapes 
are estimated, and this may also lead to a bias in the 
measured slopes. Indeed, the complexity of the exercise 
is illustrated by the significant differences in a m some- 
times inferred by different authors for the same galaxies, 
even when the same raw spectra are used (dB03; S03). 
To meaningfully compare the a m obtained to the cuspy 
c&int predicted by CDM, then, a thorough assessment of 
these potential biases is required. 

However, there is no clear consensus among those who 
have undertaken this task. The simulations of dB03 
indicate that no single systematic effect can reconcile 
their data with cuspy CDM halos, but those of S03 and 
Rhee et al. (2004) lead them to conclude that a range 
;$ cant ^ 1 is consistent with the a m derived. With 
the present body of data, it is therefore unclear whether 
the halo shapes measured in this manner signal a genuine 
conflict for the CDM paradigm, or whether they can be 
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Fig. 2. — Distribution of a) D, b) i, c) V ro t, and d) fi^ for the 
galaxy sample. 



reconciled by taking observing and data reduction tech- 
niques into account. In the present paper we address 
this issue by deriving a m for 165 low-mass galaxies un- 
der assumptions of spherical symmetry and a minimal 
disk, using techniques similar to those of dBOl, dB03 
and S03. We then perform simulations to determine the 
impact of the observing and data processing methods on 
our results. The size, homogeneity, and parameter space 
probed by the sample allow for both an independent mea- 
surement of the distribution of a m in low-mass systems 
and a thorough investigation of the potential biases on 
the a m values obtained. 

The organization of this paper is as follows: our sample 
selection criteria are given in § and our analysis of the 
sample RCs and the values of a rn obtained are in §[21 The 
simulations we perform to compare the distribution of 
a rn with those expected for model galaxies with various 
aint are described in § and we compare our simulation 
results to the sample data in § El A summary of our 
findings as well as a list of caveats regarding the accurate 
determination of inner halo shapes is in §6. 

2. SAMPLE SELECTION 

The dwarf galaxy sample is selected from the SFI++ 
Tully-Fisher database maintained at Cornell University, 
a compilation of ^4800 galaxies with accurate I-band 
photometry and estimates of rotational velocity V ro t from 
either single-dish H I or optical long-slit spectra. The 
SFlH — h includes datasets published by our group (Dale & 
Giovanelli 2000; Vogt et al. 2004 and references therein) 
and the southern spiral galaxy sample of Mathewson and 
collaborators (Mathewson et al. 1992; Mathewson & 
Ford 1996) with photometry reprocessed by our group, 
as well as recently obtained data. For this study we se- 
lect systems from the subset of galaxies in the SFlH — h 
with archived Ha RCs. 



To obtain a sample of galaxies for which the minimum 
disk assumption is likely valid, we first select by mass, 
requiring that each system have V ro t < 130 km s -1 . This 
yields 376 dwarf galaxy candidates with high-resolution 
Ha RCs. The adopted upper limit on V ro t is arbitrary 
in the SFI++, but roughly corresponds to the threshold 
below which disk galaxies lack prominent dust lanes (Dal- 
canton et al. 2004), have little or no bulges (Kauffmann 
et al. 2003), and fall systematically below the Tully- 
Fisher relation (Kannappan et al. 2002). We then ex- 
amine the optical and I-band emission of each of these 
candidates, and discard those with distorted morpholo- 
gies such as strong warps or other irregularities. We also 
eliminate systems with evidence for a strong bar or other 
pathologies deduced from variations in I-band ellipticity, 
surface brightness, or position angle (PA) across the disk 
(Haynes et al. 1999, hereafter H99). We select against 
the presence of bulges by requiring that the measured to- 
tal I-band central surface brightness be no more than half 
a magnitude greater than the extrapolated I-band disk 
central surface brightness. We do not correct for internal 
extinction, which for dwarf systems should be quite low 
(e.g. Giovanelli et al. 1995; Giovanelli & Haynes 2002). 
We then examine each galaxy in the Two Micron All Sky 
Survey Extended Source Catalog and eliminate systems 
with irregularities or strong bulge contributions not ev- 
ident in the I-band. These criteria reduce the number 
of low-mass candidates from 376 to 268. Finally, we dis- 
card systems for which the inner RC points do not yield 
a reliable estimate of the inner halo density because of 
asymmetries, undersampling, or other distortions, yield- 
ing a sample of 165 systems. 

In principle, one could verify the assumption of a min- 
imum disk by producing multi-component mass models 
to assess potential contributions of the disk and the bulge 
to the observed RCs. This is impractical for the present 
sample, however, because the outer halo shapes (and 
hence their contribution to the observed RCs) are largely 
unconstrained by the optical RCs alone. The available 
I-band photometry constrains the stellar distribution as 
a function of r across the optical disk. However, the 
normalization of that distribution may vary by factors 
of a few given the uncertainties in the stellar mass-to- 
light ratio (e.g Bell & de Jong 2001), and for this sample 
the degeneracy between stellar and halo contributions is 
exacerbated by the lack of outer halo constraints. Mass 
models are therefore of little use in determining the bary- 
onic contribution to the optical RC for the systems in our 
sample, and we do not construct them. 

Fig. G] shows the sky distribution of the galaxies in the 
sample, and Tabled presents their individual properties. 
The majority of the galaxies in our sample (104/165) 
have long-slit optical spectra from the southern sky sur- 
vey of Mathewson and collaborators, hereafter MFB 
(Mathewson et al. 1992; Mathewson & Ford 1996), who 
included a larger proportion of dwarf systems in their 
samples than other contributions to the SFI++. These 
RCs have also been folded and deprojected by Persic & 
Salucci (1995), and thus our sample overlaps with theirs 
by the same amount (see also Persic et al. 1996). All the 
galaxies in the sample were observed with a 2" - wide slit. 
Fig.|21shows the sample distributions of distances D, in- 
clinations i, V ro t and central I-band surface brightnesses 
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TABLE 1 

Sample Properties and Measured Inner Halo Slopes 
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2272 


34 


12.2 


29.3 


17.4 


66 
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11.5 


-0.05 


0.23 


0.11 


7687 


12 32 36.5 


+02 39 37 


1715 


34 


13.0 


63.0 


17.5 


90 


56 


53 


22.2 


0.04 


0.12 


0.05 


8067 


12 57 12.0 


-01 42 24 


2833 


50 


12.5 


44.9 


19.1 


81 
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8.3 
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0.14 


0.08 


8924 
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3591 


59 


13.6 


39.4 
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90 
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0.00 


0.29 


0.11 
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14 21 13.2 


+03 26 08 


1479 


28 


11.0 


76.9 


18.4 


80 


127 


117 


10.3 


0.01 


0.12 


0.05 


9888 


15 33 05.7 


-01 37 42 


2805 


42 


12.6 


31.0 


19.3 


52 


127 


115 


5.3 


0.02 


0.20 


0.12 


10641 


16 58 06.4 


+58 53 08 


5241 


78 


14.8 


32.5 


19.5 


90 


112 
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9.0 


0.01 


0.19 


0.14 


12641 


23 30 54.6 


+20 15 02 


2713 


39 


14.0 


31.4 


20.0 


78 


89 


68 


6.2 


0.06 


0.30 


0.07 


20471 


00 42 14.7 


-18 09 39 


1536 


27 


13.1 


70.0 


19.6 


90 


85 


77 


20.4 


0.04 


0.17 


0.08 


20879 


01 11 12.5 


-58 50 20 


4817 


73 


14.8 


30.3 


20.1 


90 


94 


94 


8.3 


0.01 


0.14 


0.09 


20893 


01 12 19.0 


-45 53 17 


7980 


115 


14.0 


21.0 


19.3 


77 


121 


85 


3.0 


0.06 


0.37 


0.22 


20948 


01 15 52.4 


-54 33 15 


5450 


82 


13.5 


25.4 


19.2 


71 


118 


108 


5.1 


0.02 


0.46 


0.21 


21179 


01 32 25.7 


-34 36 30 


5304 


80 


13.4 


20.2 


18.0 


67 


104 
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8.1 


0.03 


0.27 


0.15 


21214 


01 32 48.4 


-79 28 25 


1764 


29 


11.6 


54.1 


19.5 


60 


98 
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35.7 


0.07 


0.08 


0.03 


21281 


01 36 33.7 
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4107 


63 


14.6 


42.4 


20.2 


90 
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147 


35.7 


0.00 


0.06 


0.04 



Note. — The complete version of this table is in the electronic edition of the Journal. The printed edition contains only a sample. 
(1) Uppsala General Catalogue (Nilson 1973) number where available, or else our private Arccibo General Catalogue database number. (2) — (3) R.A. 
and Dec, from the literature or measured by us on the POSS-I. (4) Heliocentric (optical) radial velocity measured from best-fit URC curve to long-slit 
data. See §|^j (5) Distance, computed from Vq and (a, 5) using the multi-attractor flow model of Tonry ct al. 2000, with H a — 70 kms -1 Mpc -1 . 
(6) Total apparent I-band magnitude mj, extrapolated to 8 optical scale-lengths (H99). (7) Optical radius r op ±, corresponding to the radius 
encompassing 83% of the total I-band light (H99). (8) Extrapolated I-band disk central surface brightness, corrected for inclination (H99). (9) 
Inclination obtained from the I-band disk ellipticity (H99) assuming an intrinsic disk axial ratio of 0.2. (10) Rotational velocity, taken as the 
velocity of the best-fit polyex curve (eq. (2j) at r op t- (11) - (13) Parameters of best-fit the polyex curve (cq. J2j) to the RC points. (14) - (15) 
Measured value of inner slope and associated error Aa m . Sec §|3] 



[Iq, the latter extrapolated from the disk surface bright- 
ness profile and corrected for inclination (H99). The dis- 
tribution of fig in the sample lies in between the LSB 
and HSB Ursa Major cluster samples of Tully & Verhei- 
jen (1997). A comparison with previous studies deriving 
a m via the minimum disk assumption (dB03; S03) is 
not straightforward, as a conversion from the B- or ri- 
band values listed in these papers is required. Adopting 
a canonical conversion of B-I=0.5, however, most of the 
galaxies of S03 fall within the peak of the /1q distribution 
in Fig. |21i- The distribution of /ip thus resembles those 
in S03 and dB03 even though we did not select based on 
surface brightness, providing further evidence that the 
minimum disk approximation is valid for our sample. 

There is some disagreement in the literature regarding 
the range of i, D or M] a galaxy should have to be in- 
cluded in studies of inner halo shapes. For each of these 
attributes, the parameter space probed by our sample 
of 165 galaxies is larger than that in previous minimum 
disk studies: this allows for a systematic investigation of 
the impact of i, D and Mj on the value of a m returned. 
Accordingly, we construct a series of subsamples, labeled 
fc=l-8, from the 165 galaxies by imposing limits on i, D, 
Mi and other parameters. A description of each subsam- 
ple is given in the first 3 cols, of Table El 

Of particular relevance in the measurement of a m is 
the distribution of sample galaxy distances (Fig. ElO- 
Recent CDM simulations predicting the value of at n t re- 
solve dark matter halos to ~ 1% of the virial radius, 
corresponding to the inner 1-2 kpc of low-mass galax- 
ies (Navarro et al. 2003). The reliability of a measured 
inner halo slope thus depends on the number of resolu- 
tion elements in the sample RCs on this physical scale. 
However, the D at which a m measured via the minimum 



disk assumption fails to recover the inner halo shape is 
unclear, since it is muddled by the manner in which raw 
spectral profiles are smoothed, resampled and folded in 
the RC derivation (dBOl; dB03; S03). The inner regions 
of the most distant galaxies in this sample (D > 150 
Mpc) are poorly resolved, with only 1 or 2 RC points 
in the inner 2 kpc. For this reason we impose an upper 
limit on the range of D in half of the subsamples (3, 4, 
7 and 8; see Table EJ). Nonetheless, we also analyze the 
complete sample, which allows for a systematic investiga- 
tion of the dependence of a m on D (§3.1), and provides 
a consistency check for the simulations we perform (§4). 
The most stringent restriction in D is made for subsam- 
ple 4 (D < 50 Mpc), which includes 46 systems. The 
distribution of D therein most resembles that in dB03 
and S03, although the mean D for the samples in these 
studies is lower than that in subsample 4. There is no 
overlap between our sample and those in previous studies 
(dBOl; dc Blok & Bosma 2002; Marchesini et al. 2002; 
dB03; S03; Gentile et al. 2004). 

3. ROTATION CURVE ANALYSIS AND INVERSION 

As outlined in § El many of the RCs in our sample are 
mined from previously published studies; the reader is 
referred to those papers for details on the extraction of 
rotation velocities from the raw long-slit spectra. Previ- 
ously unpublished data have been reduced in the same 
manner as in Catinella (2005). 

Despite some differences in the techniques used to 
extract RCs from raw long-slit spectra among sample 
galaxies, the subsequent analysis was performed in a ho- 
mogeneous manner. This is fully described elsewhere 
(Catinella 2005); here we summarize the procedure. We 
fold each RC about the coordinates (x Q ff, deter- 
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TABLE 2 
Subsample Descriptions 



Subsamplc 


Description 


^ Galaxies 


Med fa™ ± Aa™) 


V / 


(2) 


(3) 


(4) 





all data (see § |2j 


165 


0.26 ± 0.07 


1 


i < 80° 


91 


0.27 ± 0.07 


2 


i < 70° 


39 


0.28 ± 0.06 


3 


D < 100 Mpc 


124 


0.23 ± 0.07 


1 


D < 50 Mpc 


46 


0.22 ± 0.08 


5 


Afj > -20.5 


97 


0.23 ± 0.07 


6 


\x off \<2", \V* ff \ <6kms- 1 


91 


0.27 ± 0.09 


7 


D < 100 Mpc, <5V < 7kms- 1 


73 


0.23 ± 0.07 


8 


i < 85°, D < 100 Mpc, Mj > -21, \x off \ < 3" 


46 


0.25 ± 0.07 




\V?f f \ < lOkms" 1 , SV < Skins- 1 






Note. — (1) Subsamplc 


name. (2) Subsamplc description. (3) Number of galaxies. 


(4) Median value of a m 


, and median error Aa m . 



mined via a least-squares fit of the universal rotation 
curve (URC) of Persic et al. (f996) appropriate for 
the Mi of the corresponding galaxy. The corrections 
(x f f ,V^j:) are measured relative to the kinematic cen- 
ter (xkin,Vkin) of the RC (Dale et al. 1997), where Vkm 
is the average of the velocities above which 10% and 90% 
of the RC points lie, and x^m is the corresponding po- 
sition along the slit. This URC folding procedure is su- 
perior to one about the kinematic center or about the 
peak of the continuum emission when either the RC or 
the galaxy light distribution is asymmetric (Dale et al. 
1997, Catinella 2005). The folded RCs are then least- 
squares fitted with an empirical "polyex" curve: 



V pe (r) = V (l 



-r/r p 



)(l+pr/r pe ) 



(2) 



where V sets the amplitude of the fit, r pe is a scale- length 
that governs the inner RC slope, and (3 determines the 
outer RC slope. This functional form is sufficiently mal- 
leable to model the wide variety of RC shapes found in 
the SFI++. Our sample galaxies are equally well charac- 
terized by a polyex or URC curve; we adopt the former 
merely for practical purposes, as high-quality polyex fits 
are archived for all SFI++ galaxies with Ha RCs. The 
best-fit parameters V , r pe , and (3 for each galaxy in the 
sample are given in cols. 11-13 in Tabled Throughout, 
we adopt V ro t = V pe (r op t), where r opt is the radius within 
which 83% of the I-band light is contained (H99). 

As the majority of the RCs in the sample are from 
MFB, we adjust the properties of the RCs from other 
sources to match those of the latter. Where required, 
we drop RC points so that the effective resolution of all 
the RCs is 2" before folding. As the MFB RCs have no 
formal error estimates, we adopt the mean deviation of 
the RC points from the polyex fit within a user-specified r 
as the error 8V{ri) for each RC point. For consistency, we 
use this definition of 5V(ri) even for RCs for which formal 
error estimates are available. Fig. [3] shows a selection of 
12 RCs from the sample with these uniform properties, 
chosen to illustrate the range of RC shapes and data 
quality in the sample. The best-fit polyex curves for 
each are over-plotted, V ro t is indicated by the horizontal 
line, and the AGC/UGC number is given in the bottom 
right corner of each panel. 

Under assumptions of a minimal disk and a spherical 



halo, the halo density profile can be calculated by invert- 
ing the observed RC via a direct application of Poisson's 
equation: 



p(r) = 



1 



V(r) dV(r) V 2 (r) 



(3) 



4-7T G \ r dr r 2 
We note again that it is unphysical to ignore the lumi- 
nous components completely, even in dwarf and LSB sys- 
tems (Fall & Efstathiou 1980). However, in dark-matter 
dominated systems there is very little change in a m val- 
ues when luminous components are included (de Blok & 
McGaugh 1997), and minimum disks provide upper lim- 
its on the true halo steepness (dBOl). In addition, we 
wish to explicitly test the usefulness of this technique in 
discriminating between cuspy and corelike halo shapes. 

To compute density profiles using eq. J3J), we assign 
V{ri) = V pe (ri) using the best-fit V , r pe and (3 in eq. J2J 
at each measured point in the folded RCs. Note that 
we do not resample the RCs after folding or fitting them, 
and so the may not be evenly distributed. The deriva- 
tive of the RC is also estimated from the polyex fit: 
dV(r)/dr\ ri = dV pe (r) / dr\ ri . We discard points with 
ri < 0.1" when calculating p(ri), as the fractional errors 
in V(ri) for these points are prohibitively large. 

We characterize the inner halo shape of the density 
profiles as p(r) ~ r~ Qm for small r. For each galaxy, the 
measured inner slope a m is estimated by a linear least- 
squares fit to the inner 2-3 points of a log(p(ri)) - log(ri) 
plot, and the errors input into the fit are the larger of the 
"up" and "down" errorbars on \og(p(ri)). Because of the 
large relative errors in log(p(ri)) for these inner points 
a least-squares fit is not strictly valid, and the returned 
formal errors on the fit are meaningless. We therefore 
estimate an error Aa m in a manner similar to that of 
dBOl, taking Aa m as the mean difference between the 
reported value and those obtained by including one extra 
and one fewer point in the fit. Fig. 0] shows the density 
profiles computed from the RCs in Fig. |31 We derive a m 
from the best-fit (solid) line, and points directly above 
as well as to the left of the arrow are included in the fit. 
The UGC/AGC number for each galaxy is at the bottom 
of each panel. Fig. 01 shows that in general rt r < 1 kpc, 
where rv is the outermost point included in the fit. 

3.1. Results for the Galaxy Sample 
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Fig. 3. — Selection of 12 (i-corrected) RCs to illustrate the range of sample RC shapes and data quality. The radial coordinate is 
expressed in units of the optical radius r op t, within which 83% of the I-band light is contained (H99). In each panel, the solid red line is 
a least-squares polyex fit (eq. to the data, the horizontal dotted blue line is the adopted value of V ro t, and the UGC/AGC number is 
given in the lower right corner. 



The values of a m and Aa m obtained for each galaxy 
in the sample are given in cols. 14 and 15 of Table 
and histograms of these quantities for subsamples 0, 3 
and 8 are shown in Fig. and illustrate the distribu- 
tions of a m for a variety of selection criteria. The his- 
tograms corresponding to the other subsamples are sim- 
ilar in shape, and the median values of a m and Aa m 
for each are given in col. 4 of Table [21 When all 165 
galaxies in the sample are included, the median values 
of the measured inner slope and corresponding uncer- 
tainty are ot m = 0.26 ± 0.07. The median values for the 
8 subsamples are similar, and in agreement within the 
median Aa m . The subsample with the lowest median, 
a m = 0.22 ± 0.08 (subsample 4), is that containing only 
the nearest galaxies: this is in good agreement with the 
results of dBOl and dB03, who find a distribution peaked 
at a m = 0.2 with a tail toward steeper slopes. 

For our sample, Fig. 5 illustrates that some of the 
asymmetry towards steeper slopes in the distribution of 
a m for subsample stems from the inclusion of systems 
with D > 100 Mpc, as the "end" of the tail (a m > 0.6) 
disappears in the distribution for subsample 3. There is 
no clear change in histogram morphology when further 
parameter restrictions are imposed (e.g. subsample 8). 
For all subsamples, Aa m is typically a significant fraction 
of a m ; each individual estimate of a m thus has little sig- 
nificance. The distribution of values, on the other hand, 
does place some constraints on the range of inner halo 
shapes measured. It is clear from Fig. and Table |2 
that the measured inner slopes a m are much shallower 
than the intrinsically cuspy halos ~ 1) predicted by 
CDM simulations of structure formation. The data thus 
exhibit the cusp/core problem. 

In Fig. El we show mean values of a m for the sample as 



a function of other galaxy parameters. The points show 
the bin centers with errorbars given by (jj/^/nj, where 
<Xj is the standard deviation of the mean for rij values of 
a m in the j th bin. There is a clear trend of increasing a m 
with increasing D, which we explore in more detail below. 
There is also a decrease in a m with increasing r pe : the 
steeper the inner RC rise as characterized by the polyex 
fit, the larger the measured value of a m (see § 14. 111 . In 
addition, there is some evidence that galaxies at higher 
i have lower a m , as found by Rhee et al. (2004). This 
could be the result of greater slit smearing at higher incli- 
nations, extinction effects or baryonic distortions hidden 
in edge-on disks and hence not apparent in the I-band 
photometry. There is also some indication that galaxies 
with Mi < —20.5 (as in subsample 5) have larger a m 
than fainter systems, which may signal a non negligible 
contribution from the stellar disk in the highest mass 
galaxies in the sample. However, no corresponding trend 
is found between a m and V ro t, nor between a m and /Ltg. 
Similarly, there is no obvious dependence of a rn on the 
offsets x f f or obtained during the data processing 
stage (see §|3J). 

The effect of RC resolution on a m is illustrated in 
Fig. [7| where a m ± Aa m is plotted as a function of r in , 
the location of the RC point that lies closest to r = ex- 
pressed in kpc. The solid lines show the expected slope 
a = |dlogp/dlogr| (evaluated at ri n ) for an isothermal 
model p oc [1 + (?*/f c ) 2 ] 1 often used to parameterize 
halo shapes. From left to right, models with r c = 0.5, 
1 and 2 are overplotted. The dashed lines show a ex- 
pected for NFW halo shapes (eq. that bracket the 
range of V ro t in the sample (see § 14.2.11) : from top to 
bottom, models with (c = 12, V200 = 50kms~ ), (c = 
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Fig. 4. — Density profiles for the RCs in Fig, l3l from eq. J2J. In 
each panel, the solid red line shows the linear least-squares fit used 
to derive a m . The points directly above and left of the blue arrow 
were included in the fit. The UGC/AGC number is at the bottom 
of each panel. 



8, V 200 = HOkms" 1 ), and (c = 5, T^oo = 180 kms -1 ) 
are overplotted. The dotted lines show CDM halo shapes 
(eq. ©) analogous to the middle dashed line, but with 
oiint = 1/2 (bo ttom) a int = 5/4 (middle), and a int = 
3/2 (top; see § 14.2. 2|) . Qualitatively, the distribution of 
points in Fig- El is very similar to that in fig. 3 of dBOl 
(see also fig. 14 in de Blok & Bosma 2002 and fig. 12 
of S03), despite significant differences in the manner in 
which ri n is derived. In particular, the RCs of dBOl are 
smoothed and re-sampled after folding whereas the ones 
in the present study are not. 

From Fig. [3 it is clear that the a m values in our sam- 
ple are in better agreement with the expected slopes of 
pseudo-isothermal halos than any of the NFW-type ha- 
los, and follow the expected trend of increasing a m in 
the pseudo-isothermal models for log(ri„)> —0.3. At 
small r jn , however, the measured values do not converge 
to a m ~ 0, as found by dBOl, but remain relatively 




Fig. 5. — Distributions of a) a m and b) Aa m for all sample 
galaxies (empty histogram), galaxies in subsample 3 (hatched red 
histogram) and galaxies in subsample 8 (cross-hatched blue his- 
togram). The median values of a m and Aa m for all subsamples 
are given in Table 01 



constant for —1.5 < log(rj„) < —0.3: a characteristic 
a m ~ 0.25 is measured throughout this range. The de- 
pendence of a m on ri n is clearly a resolution effect, first 
noted by dBOl. More distant systems are more poorly 
sampled than those of nearby galaxies in our long-slit 
spectra, resulting in the measurement of a m over a larger 
extent of the disk. A correlation between a m and D is 
therefore expected from the shape of p(ri) (Fig.QJ, since 
the slope of the density distribution steepens as r in- 
creases. This same phenomenon is likely responsible for 
the lower median a m ± Aa m in subsamples 3, 4 and 7, 
the tail toward larger a m for subsample in Fig. 
and the trend between D and a m in Fig. The need 
for spatial resolution in determining the halo shapes of 
galaxies is clear in Fig.0and underscores the importance 
of taking resolution effects into account when interpret- 
ing the results. The small, relatively constant value of 
a m obtained for log(rj n )< —0.3 also agrees with Salucci's 
(2001) conclusion that all spiral galaxies have large core 
radii. 

The discrepancy between the a m obtained and the 
cusps predicted by the CDM paradigm is actually greater 
than Figs. [S] and [7] imply for a number of reasons. Re- 
call that the a m were derived for minimum disks, an 
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Fig. 6. — Mean value of a m for the sample galaxies as a function 
of a) D, b) i, c) V rot , d) r pe , e) f) Mj, g) a; ^, and h) V* ff . 

The points show bin centers, and errorbars are given by aj/^nj, 
where Uj is the 
in the j th bin. 
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Fig. 7. — Inner slopes a m ± Aa m for the sample galaxies as a 
function of the innermost point ri n of the folded RCs. Points for 
which errorbars are not shown have Aa m smaller than the point 
size (see Fig. ISland Table ITt . Solid red lines show |d log p/d logr| 
evaluated at r^ n for a pseudo-isothermal halo with r c = 0. 5 (le ft 
curve), r c = 1 (middle curve), and r c = 2 (right curve; see §E3)- 
Dashed blue lines show expectations for NFW halos (eq~_Q) 
consistent with the range of V ro t in the sample (see § |4~2^1I : 
c = 12, V200 = 50kms _1 (top curve), c = 8, V200 = HOkms -1 
(middle curve) and c = 5, V200 = 180kms — 1 (bottom curve). 
Dotted lines show |d log p/d log r| for halos with similar to NFW 
(eq. JHJ), but with aj n t = 1 /2 (bottom), cti nt = 5/4 (middle) and 
a in t = 3/2 (top; see 



assumption that yields an upper limit on the true halo 
steepness. In addition, Id log p/d log r\ values plotted for 
the NFW halos in Fig. {7\ do not take the impact of the 
baryon collapse on the halo shape into account, which 
tends to increase aj n t over the values determined from 
collisionless simulations (Blumenthal et al. 1986; Florcs 
et al. 1993; Gnedin et al. 2004). Recent simulations also 
indicate that aj n t from collisionless simulations is larger 
than the NFW value for low-mass galaxies (Navarro et 
al. 2003); the dotted line corresponding to a.i n t = 5/4 
in Fig may therefore better reflect CDM predictions 
than NFW halos. Finally, a m is measured from a fit 
to points distributed between ri n and r^r ~ 1 kpc, and 
Figs. E] and [7| show that |dlogp/dlogr| > ai n t for r > 0. 
However, there is little change in the middle-dashed and 
dotted lines of Fig. for r < 1 kpc, where a m is typi- 
cally measured (Fig.0J. Comparing a m directly to cti n t 
is therefore a reasonable approximation for the majority 
of the sample galaxies, and adequate for the qualitative 
analysis performed here. 

4. SIMULATIONS 

The results of H3.ll demonstrate that in general, the 
galaxies in our sample have a m values that are much 
more core-like than the cuspy a^t values that CDM pre- 
dicts. The data are clearly consistent with a scenario in 
which low-mass galaxies have intrinsic cores rather than 
cusps. However, a variety of observational effects sys- 
tematically lower a m relative to cti n t to varying degrees 
if the latter is non-zero. Conversely, the dependence 
of a m on D evident in Figs. 03 HJ1 and [7| demonstrates 
that the properties of the sample itself may also influ- 
ence the values of a m obtained. We also do not account 
for |d log p/d log r I > a^t at th e loc ation where a m is 
measured in the comparisons of § l3.1l fFig.l?|). It is there- 
fore necessary to assess the extent of the potential biases 
and resolution effects on a m in order to understand the 
implications of our results for the CDM paradigm. 

In this section we investigate the nature of the 
cusp/core problem by simulating observations of model 
galaxies with various a^t and by comparing outputs to 
the distributions of a m obtained in each of the subsam- 
ples. These simulations have a single goal: to determine 
whether long-slit observations, in conditions typical of 
those in the SFI++ and of model galaxies with cuspy 
oti n t that are broadly consistent with the predictions of 
the CDM paradigm, may yield a distribution of a m sim- 
ilar to that obtained for the data. We require that the 
model galaxies have attributes similar to those of each 
subsample (e.g. Fig. EJ, and therefore simulations that 
recover the distribution of a m should also recover trends 
with other measurable properties (Fig. |BJ). The NFW 
profile has been extensively discussed in the literature 
(see §QJ. We therefore concentrate on the a^t = 1 case, 
both because of its relevance to other work and because 
prescriptions for generating NFW halos appropriate for 
low-mass galaxies are readily available. We then modify 
our technique to consider models with on n t = 1/2, 5/4 
and 3/2. 

We review the potential biases on a rn that we examine 
in our simulations in S I4.1I In § 14 .21 and § 14.31 we describe 
the model galaxies and mock long-slit observations used 
in the simulations, respectively. 
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4.1. Potential Biases on a m 

The uncertainties inherent in obtaining long-slit spec- 
tra along the major axes of low-mass galaxies, as well 
as the galaxy geometries themselves, may bias the shape 
of the resulting RCs. The impact that these biases then 
have on a m is illustrated by considering the relationship 
between a m and the RC scalelength r pe of the polyex 
function for our sample. From the values of (3 in Tabled 
and the form of eq. for r < r pe we can approximate: 



~ 1 - e~ r l r ^ = 1 - P(r) 



(4) 



Substituting this into equation @ and computing a — 
|dlogp/dlogr| gives: 



! [P(r) - 2P(2r 



pe 



[1 - 2P(r) + P(2r)] 



1 pe 



(P(r) - l)[x(P(r) - 1) - 2rP(r)] 



(5) 

In equation J5J), the second fraction is positive and varies 
little with r pe for r < 1 kpc. We therefore expect a m 
to decrease with increasing r pe for the sample galaxies, 
as is clearly the case in Fig. |SJi. An analogous effect 
will occur for a m measured for a single galaxy, if r pe is 
systematically increased. Eq. (0J, the shape of the RCs 
in Fig. [3] and their corresponding polyex fits (Table ^) 
show that a decrease in the amplitude and/or derivative 
of V(r) in the inner RC regions results in larger r pe . 
Qualitatively, then, a mechanism that biases V(r) or its 
derivative low relative to its true value will in turn lead 
to an underestimate of a rn . 

A number of potential RC biases have been discussed 
by S03 and dB03. Here, we review those addressed by the 
simulations described in § 14.21 and § 14.31 First, the finite 
angular resolution of the observations as well as the finite 
width of the slit lead to the inclusion of gas away from 
the major axis in the long-slit spectra. The projected 
radial velocities away from the major axis are lower than 
those along it; the gas at these locations may therefore 
lower the V(r) fit to the emission line shape along the slit. 
The same is true if the slit is not aligned with the galaxy 
major axis, with the impact on the RC shape exacerbated 
by the lower RC derivative away from the major axis. In 
addition, the RCs of highly inclined systems may have 
low V(r) and derivative because the disk geometry forces 
the inclusion of more gas at low projected velocities along 
the line of sight of the slit than for systems at lower i. 
Since all these effects bias V(r) and its derivative low 
these uncertainties will add coherently in real systems, 
and from eq. (JSJ we anticipate that a m < 0£i„t. 

The magnitudes of these effects have been analyzed 
by S03 (see their fig. 8), who measure a m for model 
galaxies with a range of distances and inclinations, and 
"observed" with a variety of slit offsets. For halos with 
&int ~ 1) Q: m depends most strongly on the degree of slit 
mis- alignment. Given that slit offsets are generally not 
directly measurable (see § I4.3|l , that each estimate of a m 
is highly uncertain (Fig.|BJ|, and that the models are nec- 
essarily simple in comparison to real systems (see S I4.2[) . 
these simulations cannot be used to reliably correct a m 
on a case-by-case basis (S03; dB03). Nonetheless, the 
impact of these biases can be assessed from a statistical 
standpoint by simulating observations in conditions sim- 
ilar to those in which the sample was derived, and using 



model galaxies with properties that reflect those of the 
sample. 

We note that the resolution effects discussed in § 13.11 
will bias a m in the opposite sense from that discussed 
above. To account for this, the distribution of D for 
the model galaxies is the same as that for the subsam- 
ples, and the same algorithms are used to derive a m for 
both the model and actual RCs. Our simulation outputs 
therefore include resolution effects in addition to poten- 
tial observing and geometric biases. 

There are also more complex factors that bias RCs 
derived from long-slit spectra but are beyond the scope 
of our simulations: examples are irregular Ha distribu- 
tions, non-circular motions, finite disk thicknesses and 
small, undetected bulges or bars in the sample systems. 
These last three effects are examined in detail by Rhee 
et al. (2004), who use N-body simulations of dark and 
stellar components to assess the difference between the 
measured V(r) and the actual circular velocity along the 
major axis. For most disk geometries, these factors also 
lower the value of a m derived relative to ct{ nt = 1. Since 
these effects are likely present in our sample but not in 
the models, our simulation outputs represent conserva- 
tive estimates of the biases inherent in measuring a m . 

4.2. Galaxy Models 

We model dwarf galaxies by embedding infinitely thin, 
uniform Ha disks in spherically symmetric halos with 
density distributions given by: 



pir) 



Po 



(cx) ai "*(l + cx) 3 ~ 



(6) 



where x — rj r2oo > c is the concentration index of the halo, 
p is a characteristic halo density and r^m is the halo 
size, corresponding to the radius within which p(r) — 
200 p cr u, with p crit the critical density for closure of the 
universe. For this family of density distributions, the 
corresponding halo circular velocity curve is: 



V h (r) = V 2 



2001 



' p,(cx) 
xp(c) 



with 



H(z) 



(7) 



(8) 



(Zhao 1996) and, by definition, V200 — ^200 • I n each 
simulation, we generate a set of halos with a variety of 
V200 and c but with the same ai nt ■ We elaborate on how 
V200 and c are chosen for model galaxies with ai nt = 1 
in § [OH and with a mt ^ 1 in § KTR 

Once the halo is specified, primary galaxy parameters 
such as D and i are chosen as random deviates of the 
corresponding distributions for the subsample to which 
the simulations will be compared. The I-band magnitude 
Mj is not explicitly required in the derivation of a m but 
provides an ansatz for the least-squares fitting routines 
in the RC folding process described in § O and as such is 
assigned to each model in the same manner. We assume 
that the presence of the disk has no impact on the halo 
shape, and do not apply adiabatic contraction corrections 
to our models. 
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4.2.1. Halos with a int = 1: NFW Profiles 

For ai n t = 1 in eq. © we recover the NFW halo profile 
in eq. QJ. In this case the halo velocity curve in eq. (7J 
has an analytic solution: 



yNFW {r) 



V* 



1 ln(l + cx) - (cx)/(l + cx) 
x ln(l + c) - c/(l + c) 



(9) 



Given the minimum disk analysis of § we assume that 
the halo is the sole contributor to the kinematics of the 
model galaxies. For each halo the input V200 is there- 
fore determined by substituting V^ FW (r opt ) = V rot into 
eq. jnj) and solving for V200 j given a relationship between 
V200 and c. 

Observationally, the characteristic halo concentration 
derived from fitting mass models to hybrid Ha+H I RCs 
of dwarf and LSB galaxies is c ~ 6 (McGaugh et al. 
2003). Since we wish to probe the relationship between 
the observed a m and theoretical CDM halo shapes, how- 
ever, we adopt a distribution of c that is consistent with 
predictions from N-body simulations. NFW show that 
the value of c for a halo depends on the mean density 
of the universe at the time of halo formation, and pro- 
vide a prescription for computing logc given the halo 
mass and a specific cosmology. We perform this com- 
putation for a cosmology consistent with recent Wilkin- 
son Microwave Anisotropy Porbe results (Spergel et al. 
2003): n m = 0.27, h = 0.7 and a 8 = 0.9. Jing (2000) 
shows that the distribution in c can be well described by 
a lognormal function for a given halo mass, and Wechsler 
et al. (2002) find a standard deviation of ai ogc = 0.12 
for present-day halos with well-constrained c. 

We therefore determine V200 and c for each model 
galaxy as follows. First, a (V ro t,r pt) pair measured for 
one of the galaxies in the subsample of interest is ran- 
domly chosen. We then solve eq. JHJ for V200, by sub- 
stituting V^ FW (r opt ) = Vrot and using the V 2 oo - logc 
relation from NFW. The value of c adopted for the halo 
is that from the relation, with a scatter in logc that is 
a Gaussian random deviate with cr\ ogc = 0.12. Typical 
distributions of input V200 and c determined in this man- 
ner are shown in Fig. 03 and the dashed lines in Fig. 
show |dlogp/dlogr| for the corresponding range of halo 
shapes. The distributions of V ro t returned by the sim- 
ulations resemble those in the data, suggesting that our 
method for selecting V200 does indeed yield appropriate 
halo shapes. In addition, changes in V200 by <50 kms -1 
affect primarily the amplitude of the observed RC (eq. El) , 
by an amount less than 25 kms" 1 at r opt ; since the inner 
halo shape should depend little on the RC amplitude, 
we expect that variations in the method adopted to de- 
termine V200 will not significantly alter our results (see 

§E3J. 

4.2.2. Halos with on n t 7^ 1 

To examine the impact of our observing and data pro- 
cessing techniques on a range of cuspy intrinsic halo 

2 Mo et al. (1998) and Seljak (2002) determine V200 - V ro t 
relations for L* galaxies by examining the observed Tully-Fisher 
relation and its scatter. The galaxies in our sample generally have 
much smaller masses than these systems, however, as well as sub- 
stantially submaximal disks, and as such these relations are not 
applicable here. 
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Fig. 8. — Typical distributions of V200 (top) and c (bottom) for 
model galaxies in simulations; see S I4.2.1l for details. 



shapes, we generate halos with inner slopes that are shal- 
lower (ai n t = 1/2) and steeper [punt =5/4 and 3/2) than 
the NFW value. For a m t = 1/2 and 3/2, eq. © has an 
analytic solution, and the halo velocity curves are given 

by 



V h 1/2 (r) 



2(ill 



1 ln(l + 2cx + 2S(cx)) - 2R(cx)/T(cx) 
x ln(l + 2c + 25(c)) -2R(c)/T{c) 

(10) 



_ 1 ln(l + 2cx + 2S(cx)) - (2cx)/S(cx) 
V 2 ~\ ~ x ln(l + 2c + 25(c)) -(2c)/5(c) ' 

(11) 

respectively, where R(x) — x(l + A/3x), S(x) — 
\/x(l + x) and T(x) — (1 + x)S(x). We also simulate 
galaxies with ai n t — 5/4, a value that is a better approx- 
imation to the inner halo slope near r*j ~ 1 kpc than the 
NFW value for low-mass halos (Navarro et al. 2003; see 
§0. For ai n t — 5/4, we solve eq. JHJ by quadrature. 

There are no prescriptions available in the litera- 
ture for determining c from V200 for galaxy halos when 
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a int 7^ 1, so we cannot use eqs. ijTUt and (fTTfl directly. 
We therefore determine these parameters by "bootstrap- 
ping" from the ai n t — 1 case under the assumption that 
the family of halos in eq. ijBJ is self-similar, such that 
the peak of the circular velocity profile (V max ) occurs 
at the same value of r (r max ) for halos with different 
ai n t ■ For each randomly drawn V rot , we assig n an N FW 
halo shape according to the prescription in § 14.2.11 this 
specifies an r max and V max . We then choose V200 an d c 
for the desired V? int (r) by requiring that its r max and 
Vmax match the specified NFW profile. Then, we scale 
all the Vfo int (r) values in a simulation by the same small 
factor (< 20%), so that the centroid of the distribution 
of V ro t measured for the simulated RCs is comparable 
to that for each subsample. For a range of ai„t near 1 



(1/2 



< 



Otint 



< 



3/2) this method generates halos with 



a variety of shapes that are broadly consistent with the 
current CDM paradigm and that have the same statis- 
tical properties as the data. The dotted lines in Fig. 
show |dlog/o/dlogr| for halo shapes consistent with the 
distribution of V ro t in the subsamples. 

We note that the RC shapes for self-similar halos with 
oz-int % 0.5 or aint ^ 1.5 are quite different from those 
with a^t — 1, and the scaling of V^ int (r) required to 
match the distribution of V ro t in the sample becomes 
prohibitively large. We therefore do not generate halos 
with a^t — 0, since the goal of the simulations is to 
compare our data with "observed" halos consistent with 
the CDM paradigm (which a^t — halos are not), and 
because such halos would have to be generated via a dif- 
ferent algorithm from the one presented here. 

4.3. Long-Slit "Observations" and Analysis 

We obtain RCs for each of the model galaxies by "ob- 
serving" them with a long slit. We adopt slit dimensions 
similar to those used to obtain long-slit spectra for the 
sample galaxies: the slit is 120" long and 2" wide, with 
0.66" pixels. The instrumental response is taken to be 
Gaussian, with a FWHM of 30kms _1 . We assume a 
seeing FWHM of 1" throughout. To derive the spectral 
profile for a model galaxy we center the mock slit on the 
inclined disk, but with some error in the slit position on 
the sky (r ff ) relative to the galaxy center, error in the 
PA relative to the galaxy major axis (APA), and error 
in the spectral center of the slit relative to the galaxy's 
recessional velocity (V ff). 

In order for our simulation to reproduce any biases in 
a m that may be present in our data, the values of r Q ff , 
APA and V a ff should reflect the observing conditions 
under which the sample RCs were obtained. However, 
direct measurements of these quantities are not available. 
This is particularly the case for r Q ff, as it depends on 
(i) the pointing accuracy of the telescope, and (ii) the 
accuracy with which the true RC center can be placed 
in the slit. The former has a well-determined value, but 
the latter will depend on a variety of factors such as the 
characteristics of the galaxy light and dust distribution, 
the band in which the photometric center of the galaxy is 
estimated prior to the observations, the disk orientation 
and the seeing. Factor (ii) is of particular importance for 
our sample, as the selection criteria in §0 generally favor 
galaxies with photometric centers that are more poorly 
estimated than those of L* spiral galaxies. 

To estimate r Q ff , we use the distribution of offsets x tf 
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Fig. 9. — Variation of a) |i //| and b) |V^.j.| with D for all sam- 
ple galaxies. Small blue dots show individual datapoints, and large, 
connected red circles show the median value at the corresponding 
bin center. Errorbars on the large points show the interquartile 
range. By definition, all points must lie above the dashed line. 
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10. — Distributions of a) x ff, b) Voff anc ^ c ) ^ ^ or sam Pl e 
RCs (dashed lines) and the corresponding simulated RCs (solid red 
lines), d) APA(i) (solid line) with uncertainty apA_{i) (dotted line) 
as a function of i. In the simulation, APA(i) is a Gaussian random 
deviate with mean APA(i) and standard deviation crpyi(i). 



and Vgff for the galaxy sample, determined via a best-fit 
URC curve during the data processing stage. As outlined 
in § |3 x a f f and V^tt are measured relative to the galaxy 
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kinematic center, which itself results from the average 
of the velocities above which 90% and 10% of the RC 
points lie. The parameters x ff and V^f* are thus mea- 
sures of RC asymmetry. In real systems this asymmetry 
may stem from irregularities in the RC itself (due to a 
strong bar, patchy Ha emission, or differing RC extents 
on either side of the nucleus, for example), as well as 
from a misplacement of the true RC center and galaxy 
major axis relative to the slit. Since we select against 
the presence of RC distortions in the construction of the 
dwarf galaxy sample, their impact on x ff and V^t* is 
secondary to the effects of slit alignment errors. We note 
that the values of these offsets reflect both the point- 
ing accuracy of the telescope and the ability to place 
the RC center in the slit for the reasons stated above. 
The relation between a single slit alignment error and 
the resulting x ff and V^tf ^ s nontrivial, as it depends 
on the direction of the galaxy displacement relative to 
the galaxy major axis and to the long axis of the slit 3 . 
For the sample as a whole, however, we expect the dis- 
tributions of x ff and V^j to reflect the value of r ff 
corresponding to the long-slit observations. 

We also look for trends between |x //| 5 ano - * ne 

galaxy properties used in constructing the subsamples of 
TableEl As an example, Fig.Elplots \x ff \ and \Voff \ f° r 
all sample galaxies as a function of D, a property that 
correlates well with a m (Figs. Ek and 01 . There is some 
evidence for a larger median |x //| in nearby systems, 
likely because of their larger size relative to the slit. The 
trend is weak, however, and is not evident in |V^|. We 
find no correlation between x a f j , V^s an d other galaxy 
properties, or in the value of a m obtained (see Figs. Ek 
and Eli). We therefore use the distributions of x jf and 
V„ff for all sample galaxies to select the simulation in- 
puts, noting that this results in a conservative estimate 
of r ff for samples restricted to nearby galaxies. We use 
the same input parameters for subsamples 6 and 8 as 
well, as only the outlying |x //| and |V^^| are pruned 
therein. 

We choose APA, r ff and V ff as follows. For each 
galaxy, we conservatively assign APA according to the 
i-dependent PA error relation derived from ellipse fits to 
the SFI++ I-band data (H99) shown in Fig. HHH . such 
that APA(i) is a Gaussian random deviate with mean 
APA(i) and standard deviation 07^(1). Here APA(i) 
may either increase or decrease the slit PA relative to 
the galaxy major axis. We assume that r Q ff is in a ran- 
dom direction on the sky, and choose it as a Gaussian 
random deviate with standard deviation a ro ff. We se- 
lect <j ro f f so that the distributions of x f / and V*ff re- 
covered during the analysis of the simulated observations 
resemble those obtained for the sample data, and adopt 
°~roff — 1-3" to this end. Fig. ^| compares typical dis- 
tributions of x f f and V^j for the sample data and the 
corresponding simulated observations with this choice. 
The value of V f / should have no effect on the RCs them- 
selves, as a constant velocity offset is subtracted from the 

3 In particular, we note that a displacement along the length of 
the slit will not change x j j and irrespective of APA. These 

offsets thus do not reflect the absolute position of the galaxy in slit 
coordinates. 



raw spectrum when the kinematic center is determined. 
For concreteness, we choose V ff as a Gaussian random 
deviate with standard deviation o~ vo ff — 5 krns" 1 . 

To compute the raw spectrum, we define {x, y) as slit 
coordinates with x along the long axis and (x',y') as 
galaxy coordinates with x' along the major axis. The 
distribution of measured velocities in each pixel of the 
spectrograph is then obtained by summing the parts of 
the projected galaxy disk that fall within the slit: 

P(x, V) = S x S a R(r opt )S[V - ViiCas', y% (12) 

where H(r opt ) = l for |r| < r opt and otherwise, 5 is the 
Dirac-delta function, and Vtt(x',y') is the component of 
V£" n * (x r , y') along the line of sight. The sum in x is 
over half a seeing disk and that in y is over the width 
of the slit. The grid spacing in slit coordinates is 0."1 
in all simulations. We then construct the spectral profile 
by smoothing with a Hanning window and taking every 
third pixel as an RC point, yielding an effective resolution 
of 2" . A turbulent term is added to each point such 
that the distribution of SV in the simulations matches 
that for the data (Fig. 110b). At this stage each profile 
is examined, and only those of quality comparable to or 
better than the sample RCs are further processed. For 
the subsample distributions of V ro t, SV, D, and i from 
which the galaxy parameters are chosen, the simulated 
RCs are further processed ~ 90% of the time. These 
RCs are then folded and a m for each simulated galaxy is 
derived using the same algorithms as described in § El 

5. COMPARING SIMULATIONS AND DATA 

The large number of galaxies in our sample allows for 
a quantitative comparison of a m obtained for the data 
with those expected from long-slit observations for a pop- 
ulation with a specified o.i n t ■ As each simulated spectral 
profile must be examined during the data analysis stage, 
there is a trade-off between the number of galaxies one 
simulates for comparison with sample data and the re- 
liability with which that comparison can be made. We 
find that a simulated population twice the size of the 
subsample in question is optimal in this case, and adopt 
this prescription throughout. 

To assess the level of agreement between a given sub- 
sample and a corresponding simulation, we compute a 
Kolmogorov-Smirnov (K-S) statistic. The K-S test eval- 
uates the probability that the distributions of a m for the 
sample and simulated galaxies are selected from the same 
parent distribution. We denote the P-value of a K-S test 
between the distribution of a m for a subsample k (where 
k = — 8) and that output from a simulation with a 
given ai n t as P(k,a int ). We find that P(k,a int ) may 
vary by factors of ~ 2 when the simulations for a given 
subsample are repeated. We compare our data with sim- 
ulation outputs for ai n t = 1 in § 15.11 and for simulated 
halos with different values of cti nt in § 15.21 

5.1. Simulations with ai n t — 1 

We simulate long-slit observations of a series of dwarf 
populations, with dint = 1 and observed properties that 
match those in each of the 8 subsamples defined in § El 
For each subsample k, P(k, 1) is listed in col. 3 of Ta- 
ble El Fig. ^2 compares the distributions of a m obtained 
for subsamples (all galaxies) and 3 (D < 100 Mpc) 
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TABLE 3 



K-S Test Results for Simulations with Different a int 


Subsample 


P(k, 1/2) 


P(k,l) 


P(fc,5/4) 


P(fc,3/2) 


(1) 


(2) 


(3) 


(4) 


(5) 





0.05 


0.4 


0.001 


1E-16 


1 


0.2 


0.2 


0.2 


6E-15 


2 


0.08 


0.5 


0.002 


5E-12 


3 


0.3 


0.07 


0.006 


3E-13 


1 


0.04 


0.2 


0.3 


3E-5 


5 


0.7 


0.08 


0.05 


2E-10 


6 


0.08 


0.9 


0.6 


1E-6 


7 


0.05 


0.3 


0.1 


4E-8 


8 


0.005 


0.9 


0.4 


1E-5 



NOTE. — (1) Subsample name. (2) - (5) P-value P(k, aint) of a K-S 
test between subsample k and a simulation with inner slope aint. 



with those obtained for simulated populations with the 
same statistical properties as each of these subsamples. 
Qualitatively, there is reasonable agreement between the 
simulated and actual a m for both subsamples. The sim- 
ulations also recover the same trends in a m with simu- 
lated galaxy properties as found for the sample galaxies 
(see Fig. |BJ: we show the results for all sample galaxies 
in Fig. El to demonstrate the correspondence between 
the simulations and the data over the broadest range of 
parameter space. 

In Fig. ^JJ the subsample containing all the galaxies 
(subsample 0) is better recovered than that containing 
only the galaxies with D < 100 Mpc; the latter sim- 
ulation produces a longer tail towards large a m than 
do the data, and underpredicts the number of galaxies 
with small a m . Quantitatively, the K-S statistics com- 
puted from the cumulative probability distributions inset 
in Fig. [TT] differ substantially, with P(3, 1) = 0.07 and 
P(0, 1) = 0.4. The pairs of distributions in Fig. 1771 thus 
provide an indication of the sensitivity the K-S test to 
differences between input distributions of a m given the 
sample characteristics. 

There is no clear correlation between the selection cri- 
teria for each subsample k and the value P(k, 1) returned 
in Table |3 Subsample 4, for example, has both the low- 
est median a m ± Aa m and the most stringent D cut- 
off, but P(4, 1) lies in between those of the distributions 
shown in Fig. 1111 A direct comparison between P(k, 1) 
for the different subsamples is somewhat muddled, how- 
ever, by the different numbers of galaxies therein: in 
particular, the subsamples with the strictest selection 
criteria are the smallest, and thus more extreme dif- 
ferences between the actual and simulated distributions 
of a m are required to produce a low P(fc, than in 
larger subsamples. Indeed, the two subsamples with the 
lowest P(k, 1) are subsamples 3 (D < 100 Mpc) and 5 
(Mi < —20.5), which have the largest number of galaxies 
other than subsample 0. 

Even for the larger subsamples with P(fc, 1) at the low 
end of the range 0.07 < P(k, 1) < 0.9, there is little 
evidence for a statistically significant difference between 
the distributions of a m for the sample and correspond- 
ing simulations. In other words, the a m values obtained 
for the sample RCs are consistent with a scenario in 
which low-mass galaxies have intrinsically cuspy halos 



with ai n t = 1- We note, however, that our data do 
not require that low-mass galaxies have cusps instead of 
cores. 

It is important to verify that the agreement between 
the distribution of a m output from the simulations with 
that for the data is robust to changes in our measure- 
ment technique. To illustrate this, we note that the 
values within a seeing disk of r = are highly uncertain 
because of potential RC folding errors, in addition to the 
large relative errors in V(ri) noted earlier. Accordingly, 
we re-measure a m for the galaxies in subsample 3, this 
time discarding all RC points with < 1.5" . We then 
simulate a galaxy population and analyze the spectral 
profiles obtained with this same adjustment to our tech- 
nique. The results of this exercise are shown in Fig. 1131 
The median values of a m for the simulations and the 
data are larger than in Fig. Illb : the method by which 
a m is measured clearly affects the values obtained. As 
before, however, the distribution of simulated a m resem- 
bles that for the data. Similarly, we verify that changes 
in the median input V200 of 50 kms -1 and corresponding 
adjustments to c do not significantly alter the distribu- 
tion of a m output from the simulations. 

In summary, we can recover the shallow distribution 
of a m found in the data by simulating long-slit observa- 
tions of cuspy (ai n t = 1) dwarf halos. We find a range 
0.07 < P(k, 1) < 0.9 among the subsamples, with the 
largest corresponding to subsample and the smallest 
to subsample 3. The simulation results are robust to 
changes in the data analysis technique and to small vari- 
ations in the halo shapes input in the simulation. We 
thus find that while individual galaxies show strong de- 
partures from ai n t = 1> the sample a m are statistically 
consistent with a^t — 1 halos, when observing and data 
processing biases are taken into account. 

5.2. Simulations with ai„t 1 

For each subsample of galaxies, we carry out a se- 
ries of simulations for which a^t 7^ 1- The main re- 
sults are given in Table [21 where we list P(k, 1/2) in 
col. 2, P(fc,5/4) in col. 4, and P(fe,3/2) in col. 5. 
In Fig. El we plot the distributions of a m for all 165 
galaxies (subsample 0; left column) as well as for galaxies 
with D < 100 Mpc (subsample 3; right column) and the 
corresponding simulation outputs for populations with 
(%int = 1/2 (panels a and b), 5/4 (panels c and d), and 
3/2 (panels e and f). Qualitatively, the distribution of 
a m for simulated halos with a^t = 1/2 or 5/4 for both 
subsamples and 3 are again similar to the distribu- 
tions found in the data. The histograms corresponding 
to a-int = 1/2 are slightly shallower than that for the 
data, however, and for ai n t = 5/4 there is a larger frac- 
tion of simulated systems in the tail extending to a rn ~ 1. 
There is also a larger variation in P(k, 1/2) and P(k, 5/4) 
among subsamples than for the a^t = 1 case, with 
0.005 < P(fc;l/2) < 0.7 and 0.001 < P(fc;5/4) < 0.6. 

In general, there is a good correspondence between 
a m for the simulated populations and for the subsam- 
ple data when a^t — 1/2: the larger subsamples have 
P(fc,l/2) > 0.1 and all but one have P(fc,l/2) > 0.01. 
In particular, the values of P(k, 1/2) for subsamples 3 
and 5 are larger than the corresponding P(fc,l). This 
is not the case for all subsamples with low median a m , 
however: both subsamples 4 and 7 have lower P(k, 1/2) 
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Fig. 11. — Comparison between actual (dashed lines) and simulated (solid red lines) a m for a) all sample galaxies, and b) subsample 
3, where D < 100 Mpc. In both a) and b), the main panel shows the distributions of a m , and the corresponding cumulative probability 
distributions are inset. The P-value P(k, 1) of a K-S test is also indicated. 



0.5 

i °- 4 
0.3 

0.2 



-(a) 




-(b) 




I , i i , | i i- 


: k /J 


— ^ 


] 


[ r Y 










1 , i , i 1 , ,: 



20 40 60 80 100 120 
D (Mpc) 



1 (deg) 



0.5 


L(c) 




1 1 1- 


6 0.4 








e 








0.3 


:, , , 






0.2 




i , ,E 




80 


100 
v r„i (km s-') 


120 


0.5 


L(e) 






0.4 

E 








a 

0.3 
0.2 










-2 - 


-1 1 


2 




0.5 - 

1 0.4 - 

0.3 - 

0.2 r 



(0 



x oft (asec) 



-10 -5 5 
V*„ (km/s) 



Fig. 12. — Mean value of a m for a simulated population of halos 
with a-i n t = 1 as a function of various galaxy properties (solid red 
lines), superimposed on the corresponding relations for subsample 
(dashed lines). Plot details are as in Fig. Irjl 



than P(fc,l). Also in contrast to the simulation results 
for otint — 1, the smallest subsamples (2, 4 and 8) are 
the ones with the most extreme differences between the 
simulated and actual distributions of a m . 

The simulations with ct!j nt = 5/4 tend not to recover 
the observed distribution of a m as well as either the 
otint — 1/2 or oii n t = 1 cases. This is evident in the cu- 



mulative distributions inset in Figs. 114b and d, in which 
the simulations for both subsamples and 3 underpre- 
dict the fraction of halos with 0.3 < a m < 0.5. While 
P(k, 5/4) > 0.1 for over half of the subsamples, the latter 
tend to be the smallest ones considered, and thus more 
extreme differences between the simulated and actual a m 
are required to lower P(k, 5/4) than for the larger sam- 
ples. There is no clear trend between the median a mi se- 
lection criteria or size of a given subsample and P(k, 5/4), 
however. We thus find that long-slit observations of halos 
with OLi n t = 5/4 yield distributions of a m that are only 
marginally consistent with the sample galaxies. Sample 
sizes larger than those examined here are needed to fur- 
ther investigate this case. 

For otint = 3/2, our mock long-slit observations of 
galaxy models (Figs. 114b and d) yield a distribution of 
a m that is clearly discrepant from that found in the data. 
We do find some small a m despite the steep input o-inu 
however, most of the values are too large to be consis- 
tent with the core-like distributions in the data. As such, 
P(k, 3/2) << 1 for all subsamples, and the hypothesis 
that the distributions of a m for the data and oti n t =3/2 
simulations stem from the same parent distribution is 
strongly rejected. 

6. DISCUSSION AND SUMMARY 

Under the assumptions of a minimum disk and spher- 
ical symmetry, we have measured the inner halo shapes 
a m for a sample of 165 dwarf galaxies from long-slit opti- 
cal spectra. We find a median value a m — 0.26±0.07 for 
all sample RCs, and similar values for various subsamples 
of the data (Table EJ) . Irrespective of the subsample in 
question, the distribution of a m is significantly shallower 
than the cusps (ai n t) predicted by the CDM paradigm 
(e.g. Figs. El EJ). Our results resemble those obtained in 
other studies in which similar analyses were performed, 
albeit with smaller samples (e.g. dBOl; dB03; S03): it 
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Fig. 13. — Same as in Fig. Illi t. but with all points rj < 1.5' 
excluded in both the sample and simulation analyses. 



thus seems that in general, halo shapes inferred via the 
direct inversion of long-slit spectra (eq. ©) exhibit the 
cusp/core problem. The variation of a m with n n (Fig. 01 
clearly demonstrates the need for high-resolution obser- 
vations in determining a m from RC shapes. 

To elucidate the relationship between a m and vari- 
ous intrinsic inner slopes cti n % that are consistent with 
the CDM paradigm, we simulate long-slit observations of 
model galaxies with the same global properties as each 
of our subsamples. The clear dependence of a m on D 
(Figs. Ek, 13 as well as on the measurement technique 
fFig. I13fl . and the wide range of a m measured for model 
galaxies in these simulations (Fig. I14J) underscore the im- 
portance of simulations in assessing the nature of the 
discrepancy between a m and . 

Our simulations show that populations of dwarf galax- 
ies with a int — 1/2 or 1 observed with a long slit yield 
distributions of a m that are broadly consistent with those 
obtained for the sample RCs (Figs. fTTl I14b -b). These 
simulations also recover trends in a m with other galaxy 
properties fFig. ITU), suggesting that they capture the es- 
sential processes that govern the value of inner halo slope 
measured. The distribution of a m in our sample is there- 
fore consistent with halo profiles that are shallower than 
or comparable to the NFW case: the discrepancy be- 
tween a m and a.;„t ~ 1 in our sample is reconciled when 
the impacts of observing and data processing techniques 
are considered. 

We emphasize that this result does not imply that the 
halos in our sample actually have cusps with ai n t ~ 1: 
the available data are merely statistically consistent with 
them. Indeed, it has been repeatedly demonstrated (e.g. 
dB03; S03) that corelike distributions of measured in- 
ner slopes such as those found here are also consistent 
with a variety of intrinsically corelike halo shapes. More- 
over, our findings clearly do not address the nature of the 



cusp/core problem obtained from other types of analy- 
ses (see §0. However, they do demonstrate that, for a 
large sample of dwarf galaxies, the distinction between 
intrinsic ai n t = 1 cusps and cores cannot be made by 
measuring a m under assumptions of minimal disks and 
spherical symmetry 

We find that simulated observations of model galax- 
ies with ai n t = 5/4 are only marginally consistent with 
the data in our sample fFig. iTlb-dl. It is therefore un- 
clear whether observing and data processing techniques 
are sufficient to reconcile a m with this case, and the 
cusp/core problem may remain should CDM simulations 
converge to halos with aint ~ 1.3 on galactic scales (e.g. 
Navarro et al. 2003). A sample size larger than the 
one considered here may offer firmer conclusions for the 
a int =5/4 case. More sophisticated simulations that in- 
clude the effects discussed by Rhee et al. (2004) should 
also be considered in this case. The hypothesis that 
dwarf galaxies have aint = 3/2 is strongly rejected by 
our data (Figs.H"4*k-f). and we do not expect that adding 
complexity to our models will change this result. 

For dark matter-dominated systems such as the ones 
considered here, there is typically little change in a m 
when the luminous components are included in the anal- 
ysis (dBOl, dB03, S03). In addition, our simulations 
reproduce trends in a m with various galaxy properties 
found in the data (Fig. 112(1 . indicating that they do 
capture the processes that govern the value of a m ob- 
tained. For these reasons, we expect that the assump- 
tion of a minimum disk does not significantly affect our 
results. Indeed, the extra parameters introduced in mass 
model fits by including luminous components renders the 
O-int = 1 case prohibitively difficult to rule out, even in 
the absence of the RC uncertainties addressed here (S03). 
It is therefore unclear whether the inner regions of long- 
slit Ha RCs should be used to address the cusp/core 
problem at all: studies that exploit the RC shape near 
r opt rather than at r ~ should provide more robust 
constraints on the halo profile (e.g. Salucci 2001). 

High-resolution two-dimensional analyses of spiral 
galaxy velocity fields are immune to many uncertainties 
that plague those from long-slit spectra, and may there- 
fore be able to distinguish various ai n t from the mea- 
sured a m . These observations are underway: Bolatto et 
al. (2004) and Blais-Ouellette et al. (2004) present RCs 
for a total of 11 systems derived from Ha + synthesis 
CO maps and Fabry-Perot Ha maps, respectively. The 
5 galaxies of Bolatto et al. have different halo struc- 
tures ranging from 0.3 < a m < 1.1 The data from Blais- 
Ouellette et al. are thoroughly analyzed by Dutton et 
al. (2003), who run an extensive set of mass models with 
contributions from stellar and gaseous disks and allow for 
halo shapes considerably more complex than those con- 
sidered here. They find that ai n t — provides a better 
fit to the data than halos with ai n t = 1 for a range of 
model parameters. 

provides the best fit for all but one of their 6 RCs, but 
cannot conclusively rule out ai„t > 1 systems. 

As the quality of spiral galaxy kinematics derived from 
observations increases, the presence of baryons in dark 
matter halos cannot be ignored in CDM paradigm tests. 
The discrepancy between cusp and core in mass models 
increases if the disk contracts adiabatically in the halo 
(dBOl; Dutton et al. 2003; S03). However, the implicit 
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Fig. 14. — Comparison between distributions of a m for subsample (all galaxies; left column) and for subsample 3 (D < 100 Mpc; right 
panels) and correspond ing simulation outputs for eei nt = 1/2 (a, b), Oj n t = 5/4 (c, d), and Oj n t = 3/2 (e, f). In each panel, the plot details 
are the same as in Fig. 111! 



assumption that the initial angular momentum distribu- 
tions of the baryons and the dark matter are the same 
has been called into question (e.g. Bullock et al. 2001; 
van den Bosch et al. 2001), and Gnedin et al. (2004) find 
that adiabatic contraction over-estimates the true degree 
of halo contraction. More dramatic gasdynamical effects, 
such as star formation and feedback during disk forma- 
tion, also need to be considered. The impact of central 
starbursts during disk formation was simulated by Mo & 
Mao (2004), who find lower c and a,„( relative to those 
of collisionless halo shapes. As such, it is perhaps not 
surprising that the halo shapes inferred in regions with 
strong baryonic contributions to the mass density do not 
match the profile shapes of the collisionless, precollapse 
halos discussed here (Weiner et al. 2001, Binney 2004). 
Indeed, the discrepancy between CDM theory and ob- 
servations is largest at low redshift, where gasdynamical 
effects are most important: the role of these effects in 
galaxies and cluster - sized halos needs to be assessed 
before a CDM crisis is advocated. 
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